Electrical quantum measurement of a two level system at arbitrary voltage and 

temperature 
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We calculate the noise spectrum of the output signal of a quantum detector during continuous 
measurement of a two- level system (qubit). We generalize the previous results obtained for the 
regime of high voltages (when eV is much larger than the qubit 's energy level splitting A) to the 
case of arbitrary voltages and temperatures. When V ~ A the output spectrum is essentially 
asymmetric in frequency, i.e., the output signal is no longer classical. In the emission (negative 
frequency) part of the spectrum the peak due to the qubit's coherent oscillations can be 8 times 
higher that the background pedestal. For V < A and T = the coherent peaks do not appear at 
all. 
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I. INTRODUCTION 

The problem of quantum measurements has been 
around since the early days of quantum mechanics. The 
recent upsurge in the interest to the quantum comput- 
ing made it necessary to investigate the properties of the 
real physical systems used as quantum detectors. Thus 
far the most developed and successful mesoscopic detec- 
tors are the under-damped SQUIDs (or current biased 
Josephson junctions) which perform switching (thresh- 
old, latching) measurements. These systems were inves- 
tigated for many years. Recently such measurements 
were used to confirm the coherent dynamics and ma- 
nipulations of the superconducting qubitsi* 2 *^. Other 
measuring devices operate in a smoother mode similar 
to a (linear) amplifier. Especially interesting are electro- 
meters whose conductance depends on the charge state 
of a nearby qubit. Two families of electro- meters work- 
ing in this regime have mostly been discussed in the 
literature. These are the single electron (Cooper pair) 
transistors (SET^i^i?^ 1 ?! 11 ! 12 and the quantum point 
contacts (QPCy-VV-WTi 18 ! 1 WV 2 , Some of these 
schemes have been implemented experimentally and used 
for quantum measuremento 23 i 24 i 25 i 26 ' 27 i 28 . 

In general one distinguishes between strong and weak 
quantum measurements. In the strong (von Neumann) 
measurement the meter discriminates between the states 
of a qubit on a time scale much shorter than the other 
time scales of the system, e.g., the inverse level splitting 
of the qubit. Thus the qubit's internal dynamics is irrel- 
evant and the meter determines the basis in which the 
measurement takes place (the eigenbasis of the measured 
observable). Such measurements strongly resemble pro- 
jection of the qubit's state even though everything can be 
described by treating the coupled system of the qubit and 
the meter quantum mechanically. In the opposite, weak 
measurement limit, the measurement time is relatively 
long and the state of the qubit may in principle change 
during the measurement. Yet, there is a special regime 
in which the measurements resemble projection. This is 
the so called quantum non-demolition (QND) limit. For 
example, if the coupling is longitudinal, i.e., if a spin 



(qubit) is placed in a magnetic field along the z-axis, 
and the detector is measuring the a z observable of the 
spin, and there is no extra environment capable of flip- 
ping the spin, then <r z is conserved (non-demolition) and 
can be measured even if it takes a long time. In all other 
weak coupling regimes the meter cannot extract precise 
information about the initial state of the qubit. One can 
only talk about continuous monitoring of the qubit by 
the meter in which they influence each other. Studying 
such continuous monitoring, e.g., in the stationary state, 
is still useful, as one can extract physical characteris- 
tics of the meter and of the qubit and later use them for 
manipulations, projection-like measurements or quantum 
feedback control of the qubits2£. 

Continuous monitoring in the regime close to QND was 
considered in Refill. The small deviation from the QND 
limit causes rare spin flips. Thus the current (output sig- 
nal) in the meter shows the "telegraph noise" behavior. 
In the noise spectrum of the current this translates into 
a Lorentzian peak around zero frequency. Recently the 
QND measurements have been considered in the rotating 
frame of a spin (qubit The regime far from QND was 
the main focus of RefsiiS*^. This regime is realized, e.g., 
in the case of the transverse coupling between the qubit 
and the meter (magnetic field along the z-axis while a x 
is being measured). The possibility to "observe" the co- 
herent oscillations of a qubit was analyzed. While, due 
to the unavoidable noise, the oscillations can not be seen 
directly in the output current, the spectral density of the 
current noise has a peak at the frequency of the oscil- 
lations (Larmor frequency, level splitting of the qubit). 
The laws of quantum mechanics limit the possible height 
of the peak. In the case of a 100% efficient (quantum 
limited) detector the peak can be only 4 times higher 
than the background noise pedestal^. Inefficiency of the 
detector reduces the height of the peak further. Such in- 
efficiency implies^Si 2 ^ that when the detector is used in 
the QND regime, the measurement time, i.e., the time 
needed to discriminate between the states of the qubit is 
longer than the lowest possible limit for this time, i.e., 
the dephasing time. This, in turn, means that the meter 
produces more noise than it is necessary for the measure- 
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ment, or, in other words, that some information obtained 
by the meter is not transferred to the output signal. 

All the results described above were obtained in the 
limit when the voltage applied to the measuring de- 
vice is much higher than the qubit's energy level split- 
ting eV ^> A. In particular, in this regime, the "tele- 
graph noise" peak around ui = is absent in the case 
of the purely transverse coupling, while in the intermedi- 
ate regime (between longitudinal and transverse) the two 
peaks coexist. The output noise spectrum, in the regime 
eV 3> A is almost symmetric (classical) at frequencies 
of order and smaller than A. In this paper we relax the 
condition eV 3> A. We calculate the non-symmetrized 
current-current correlator in the case of the purely trans- 
verse coupling between the qubit and the meter for arbi- 
trary voltage and temperature. At low voltages, eV ~ A, 
the output noise is essentially asymmetric, i.e., the out- 
put signal is quantum. In other words, we have to dif- 
ferentiate between the absorption (to > 0) and the emis- 
sion (u < 0) spectra of the detector (see, e.g. RefiS). 
Thus, the detector ceases to be a device able to trans- 
late quantum information into the classical one, and the 
way the output is further measured becomes important. 
The qubit produces two symmetrically placed peaks in 
the output noise spectrum. The ratio of the peaks' hight 
to the hight of the pedestal can reach 8 for the negative 
frequency peak. We also obtain a small peak at lu = 
(even for purely transverse coupling). 

The rest of the paper is organized as follows. In the 
next section we define the physical system and outline 
the computational scheme. In Section IIIII we derive 
general expressions for average current and fluctuations 
spectrum, valid for any voltage and temperature. In sec- 
tion llVl we provide results for the specific model of a qubit 
transversely coupled to a quantum point contact. In the 
appendices we establish connections and clarify distinc- 
tions between the technique developed here and the other 
existing approaches. In Appendix lAl we analyze the va- 
lidity of the MacDonald's formula for the current noise 
in our case. In Appendix iBl we analyze the applicability 
of the Bloch-Redfield approximation. In Appendix IU1 we 
outline the perturbation theory based on spin represen- 
tation by Majorana fermions. 



II. THE SYSTEM 

We study the quantum measurement process in which 
a quantum point contact (QPC) is used as a measuring 
device. These devices are known to serve as effective 
meters of charge (see, e.g., Refsi^&Sii&Sii). In this paper 
we consider the simplest limit of a tunnel junction when 
the transmissions of all the transport channels is much 
smaller than unity and is controlled by the quantum state 
of a qubit. This model has previously been used by many 
authorsi^iSiSS. The measuring properties of QPCs in a 
more general case of open channels have been studied, 
e.g., in Refsi 21 ' 22 . 
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FIG. 1: Dyson equation. 

The tunnel junction limit is described by the following 
Hamiltonian 

H = € ' c l c ' + X! €rC r Cr 
I r 

+H sys + «(cj Q e~ lk + h.c.) , (1) 

l,r 

where the tunneling amplitude f2 may in principle depend 
on any operator of the measured system (qubit). In this 
paper we focus on the case of transverse coupling between 
the meter and the qubit 



H sys = -(l/2)A(7 i 
Q = T + T\o x 



(2) 
(3) 



which is optimal for "observation" of the qubit's coherent 
oscillations. The transmission amplitudes To and T\ are 
assumed to be real positive (as the phase of the trans- 
mission amplitude does not matter in our case) and small 
(tunnel junction limit). Obviously, we only have to con- 
sider T\ < T . We do not, however, assume T\ <C T , 
i.e., we allow for the detectors with large (relative to the 
average output) response. We have also introduced the 
counting operator: e~ lk \m) = \m + 1), where m is the 
number of electrons that have tunneled trough the point 
contact. Using this trick we could in principle study the 
full counting statistics of the current. 

Following Refi2& we integrate out the microscopic de- 
grees of freedom (the electrons) in the left and the right 
leads and consider the time evolution of the reduced 
density matrix of the system. This density matrix is a 
function of the measured system's coordinates as well 
as of the variable m, i.e., p = p(mi,m2). Due to 
the translational invariance with respect to m it is con- 
venient to perform the Fourier transform p(k\,k2) = 
Smi m 2 p( to 1j m2)e~ lkimi+lk2m2 . In this representation 
the operators e ±lk in Eq. JTJ are diagonal. The mas- 
ter equation for the density matrix with the information 
about the number of electrons that have tunneled was 
used in Ref J*i. Here we do the same in the Fourier space 
(fci and fc 2 indexes) . We write down the Dyson equation 
for the propagator of the density matrix (see Fig.^l . Tak- 
ing the time derivative, one arrives^ at the generalized 
master equation 



-p{t) - L oP (t) = J dt' E(t - p(t') , (4) 

to 

with the zeroth order Liouvillian 

L =[l® iHn - iH <g> l] . 



(5) 
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Throughout the paper we assume % = e = 1. As we have 
integrated over the electronic degrees of freedom, the Li- 
ouvillian Lq as well as the Hamiltonian Hq operate in the 
direct product of the spin's (| f / |)) and the counting 
(|m)) Hubert spaces. As in the counting space the zeroth 
Hamiltonian is zero (the number m changes only due to 
the tunneling) we have Hq - 
orously one should write Hq 
to present the Liouvillian as a super-operator acting from 
the left on the density matrix regarded as a vector. This 
way of writing makes the analysis easier and is especially 
convenient for numerical simulations. As an example, we 
rewrite the product Lqp in the matrix form: 



Q. 



<XX + > 

<x + x> 



H sys (mathematically rig- 

H sys <g> 0) . We have chosen FIG ' 2: The first order approximation for the self-energy. 



where rj = 27T/9lPr- Fig. |3 shows the diagrams that lead 
to Eq. 0. 

Formally, Eq. can be solved by applying the 
Laplace transformation (that is, it becomes a simple sys- 
tem of linear equations for each value of s): 



L oP = — 7T 




(6) 



In the lowest non-vanishing approximation (second or- 
der in the tunneling Hamiltonian) for the self-energy we 
obtain 

£(*) = 

a%(t) n u _ kl u (t) n d 2 + a*_(t) ni u (t) ni k2 

+a+(t) n d k2 U (t) Ql ki + a _(t) !li k2 U (t) 

- a+ {t) n u ki u (t) n u _ ki - a _(t) n u _ kl u (t) n u kl 
-<(*) n d _ k2 u (t) n d 2 - a* (f) at u (t) ni k2 , 



(7) 



where H± k 



e ±iki n u , a d ±k = e ±ik2 Q d , n u = (n <g> 1) 



d — 



(1(g) !V) and U (t) 



,L t 



The superscripts u and 



d stand for the "up" and "down" Keldysh contours. We 
observe that all the matrix elements of the self-energy 
E(i) in Eq. Q are functions of Afc = k\ — k 2 only, i.e., 
they describe transitions which conserve m\ — r7i2- In 
particular they connect the diagonal elements (mi = m<i) 
with only the diagonal ones. Even though we could have 
multiplied the factors e ±jfcl / 2 in Eq. J7J) and express E 
as function of Afc, we keep these factors separately in 
this particular formula for future convenience. We also 
introduce the correlators 



a + (t) ee (X(t)XH0)) , 



(8) 



and 



a_(t) ee (Xt(i)X(O)) , (9) 
where = c£ c;. Their Fourier transforms are: 



and 



a+(w) = 77 (w + V) 



.(w) = »/(w-l0 



1 , LO + V 1 

— coth 1 — 

2 2T 2 



1 , cj-F 1 
- coth 1 — 

2 2T 2 



(10) 



(11) 



where 



p(ki,k 2 ,s) = U(s,Ak)p , 



£7(s, Afc) ee (s - L - E(Afc, s)Y 



(12) 



(13) 



and po is the density matrix at t = to- 

Further approximations are sometimes used to make 
the master equation Q Markovian. When the dissipative 
processes are slow in comparison with the unperturbed 
coherent (Hamiltonian) dynamics, the Bloch-Redfield ap- 
proximation is appropriate. Within this approximation 
one substitutes /3(f) -> e -io(t-t ')j3(i) in the RHS of 
Eq. Q). Tllis lea ds to E(s) — > E B r in Eq. (T^J, where 



J BR 



dtY,(t)e 



-L t 



(14) 



In our case the validity domain of this approximation ex- 
tends also to the regime when the dissipative rates are 
bigger than A, i.e., when the qubit is over-damped. Be- 
low we will see that, due to the smallness of To and Ti, 
the qubit can become over-damped only at high voltages 
or temperatures, i.e., when y > A or T > A. In that 
case, however, the self-energy E(t— t') decays on the time 
scale of order 1/V or 1/T and, thus, is Markovian for 
slower processes. As we are mostly interested in frequen- 
cies not much higher than A, we can still use Ebr. In this 
paper we will mostly employ the Bloch-Redfield approxi- 
mation, which gives very simple and transparent results. 
In Appendix [5] however, we will use the non-Markovian 
expression (|13|l to confirm that the non-Markovian cor- 
rections are small. 

As an example of the Bloch-Redfield approximation, 
let us consider the regime V ^> A and T = 0. Then, 
for \uj\ <C V one has a + (ui) = ij(V + cj), i.e., a + (t) = 
T)(V8(t)+i8'(t)), while a-(w) = 0, i.e., a<_(i) = 0. Then 
using the definition l|14|l we obtain 



{2e tAk n u n d - n u n u - n d n d } 



{e lAk n u [Lq, n d ] - e lAk !l d [Lq, «"] } 
[L ,fi"]-fi d [La,!)*]} 



(15) 
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FIG. 3: Diagrammatic representation of the correlator 
(m(t)m(t'))- 



where the "value of the delta function at zero", S Ul ,(0), 
should be understood as a constant depending on the 
high energy cut-off uj c (diverging with it). Note that 
n u Q d = Q d Q u . The last term of Eq. JTSJ) can be rewrit- 
ten as irj6u c (0) [fl 2 , ...]. Thus, this is a renormalization of 
the spin's Hamiltonian. This renormalization is usually 
disregarded either because it is small or, as in Caldeira- 
Leggett's approach, since a counter term has been al- 
ready added in the initial Hamiltonian. Substituting the 
self-energy of Eq. 115|) into Eq. 10} we arrive at the gener- 
alized master equation obtained in Ref>2 9 by other tech- 
nique. 



III. CALCULATION FOR ANY VOLTAGE AND 
TEMPERATURE 

We start by writing down the formally exact expression 
for the correlator (m(t)m(t')) for t > t' (see Fig. |3J): 

(m(t)m(t')) = Tr [m u U(t,t')m u U(t' ,t )p ] 
t t' 

+ j dhj dt 2 Tr[m u U(t,t 1 )Sm u (t 1 ,t , ,t 2 )U(t 2 ,t )p } , 

(16) 

where m w = (m® 1) = id/dki is the bare vertex rh on the 
upper Keldysh branch while 5m u (t2, t', ti) is the vertex 
correction. The importance of the vertex corrections was 
recently pointed out in Ref>i£. The Tr[...] operator is the 
trace of the density matrix (not of the super-matrices 
like £). The trace over the ki. 2 indexes is calculated as 
X-tt J*--* ^T0T ZnSih - k 2 ).... For the stationary state 
properties, which do not depend on the initial density 
matrix, the trace operator over reduces to taking 
the limit Ak — ► 0. This follows from the fact that all the 
propagators depend on Ak only. 

For the vertex correction we use the same approxima- 
tion we have employed for the self-energy J7J) (see Fig.|3}. 
One easily obtains 8m u (t\, t', t 2 ) from Eq. Q inserting 
the operator m u in all the terms between the f2 operators 
(from either side of Uq as m u commutes with Uq): 

6m u {h,t',t 2 ) 

= a* + (t 1 -t 2 )m ki U (t 1 -t 2 )m u n d k2 +... , (17) 



5m u 



l = £2 
= m 



<XX + > 

<x + x> 



FIG. 4: The first order approximation for the vertex correc- 
tion. 



where ... stand for seven more terms obtained in the same 
way from Eq. Q. We observe that 5m u (t\, t' , t 2 ) — 
5m u (t\ — i 2 ) in this approximation does not explicitly 
depend on t' . We differentiate Eq. (|16|l over t' and t 
and use the master equation p} written in the form 
dU/dt = (L + S) *U = U* (L + E) to obtain: 



(J(i)/(f )) = Tr [m u A (t _ tt) U (t ,_ to) p ] 



+Tr 
-Tr 



m u {{L + E) * U * A} (t _ tr) U (t ^t )Po 
m u {(L + £) * U} [t _ tl] {B * U} [t ,_ to) po 



(18) 



where A = <5m u — E m" and B = 6m u — m"£. The convo- 
lutions are defined as {g* f}(t-t') = L> dtig(t — ti)f(ti — 
t') and, analogously, {g* /* fi}u-t') = Li dii fl 1 dt 2 g(t — 
ti)f{tx - h)h{t 2 - i 1 ). The symbol L in Eq. JTHJ) 
should be understood as a local in time kernel, i.e., 
L (t - t') = L S(t - t' - 0). From Eqs. £7J) and 
we obtain 

A(t) = 

a+(t) n d k2 U (t) Cl! kl - a_(t) Sl d _ k2 U (t) Q kl 
-a+(t) ill U (t) Q u _ kl + a_(t) Sl u _ kl U (t) U kl , 



(19) 



and 



B(t) = 

-a* + (t)Sl! kl U (t)Sll 

-a + {t)ni x u {t)wt ki 



a*_{t) niu Q (t) n d _ k2 
a_(t) n u _ ki u (t) ni 



(20) 

Next we note the following property of the super- 
operators L , E, and A: Tr [Zq....] = 0, Tr [£....] = 0, 
and Tr [A....] = 0. This allows us to simplify Eq. [tTSjl : 

(I(t)I(t')) = Tr [[m^,A] (t _ tl) U {t ,_ to) p Q ] 
+Tr [{[m u , £] * [/ * 4} (t _ t/) C/ (i ,_ to) p 
-Tr [{[m u ,£] * E/} (i _ t r) {£ * 17} ^ p 



(21) 



and the commutators are readily calculated: 

[m u ,A] (t) = 

a+(t) < u (t) ai kl + a _(t) tti fe2 tr (t) 



(22) 
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K,%) = 

a* + (t) SVL^ U (t) fi& - at(t) Q u kl U (t) Qt k2 

a+ {t) n d k2 u (t) n* fcl - «_(*) ni fa y (t) ^ ■ 

(23) 

We observe that the super-operators A, B, [m u ,A], and 
[m u ,E] do not contain the m operators, i.e., there are 
no differentiations over left in Eq. (|21l) . Thus we 
can safely perform the limit Afc — > 0. We introduce the 
functions a(t) = a+ (t) + a- (t) and /3(f) = a+ (t) — a_ (t) 
and obtain for fci = &2 

A{t) = -/3(t)n u u (t)n u + 0(t)n d u o (t)n u , 
B(t) = -p*(t)n u u (t)n d - p(t)n u u a (t)n u , 

[m u ,A] {t) =a(t)n d U (t)n u , 

[m™,s] (t) = /3*(t)n u u (t)n d + /3{t)n d u (t)n u . 

(24) 

Finally, for the stationary state we take to — > — oo and 
obtain (for t > t') 

(I(t)I(t')) = Tr [K\%_ t , )( 3 st ] 

+Tr[{[m tl ,£]*[/*A} ( ^ t , )( S st 

-Tr[{[m",S]*[/} (t _ t0 {£*l} (oo) 

(25) 

where by definition {£? * l}^) = J dt\B(t\) = B(s = 
+0) and B(s) is the Laplace transform of B(t). The 
stationary density matrix is given by 

p st = Sf/(s, AA)/3o| a -»0,Afc->0 • (26) 

Analogously we find the expression for the average cur- 
rent: 

(I(t))=Tr[{[m«,E}*U} {t _ to) po\ , (27) 

which in the stationary regime becomes: 

(I) =Tr[{[m",E]*l} oo/ 3 8t ] . (28) 

Equations l|24l) . (|25|l . and l|28|l constitute the central re- 
sult of this chapter. They allow us to calculate the 
current-current correlator and the average current in the 
first order approximation for the self-energy and the ver- 
tex corrections. 



IV. AVERAGE CURRENT AND NOISE 
SPECTRUM 

To formulate our results in a compact way it is conve- 
nient to introduce the two following functions: 



and 



s(w) 



_ a(u>) + a(—uj) 
= 2?t 



(V + uj) , V + oj (V-w) , V-u 
- — - — - coth — — — + - — - — - coth ■ 



2T 



2T 



a(u>) 



2r) 



(V + lu) , V + UJ (V-lj) , V-w 
- coth coth ■ 



2T 



2T 



(30) 
while 



One can easily check that a{u>) 
0(cj) = rj(a(u) + V). 

For the average current in the stationary regime we 
obtain the following expression 



(I)=g V + 9l V[l- 



A a (A) 



(31) 



(29) 



where we have introduced the conductances go = T)Tq 
and gi = j]Tf. If T = the result simplifies. For V < A 
we obtain (I) = goV , i.e, no contribution of the qubit. 
For V > A we have (I) = go^+giV^l- A 2 /F 2 ). Finally, 
for V > A we obtain (7) « ff0 V + ff X V = (I u + 7j,J/2, 
where we have introduced the values of the current corre- 
sponding to the two cc-projections of the qubit: -f^n = 
?7(T ± Ti) 2 V = V(g + gi ± 2^5001). This result be- 
comes intuitively clear if one notes that the relevant fre- 
quency scale of the tunneling process is equal to V while 
the fluctuations of spin's observable a x have a charac- 
teristic frequency A. In the regime max[F, T] <§C A the 
spin is mostly in the ground state and the tunneling elec- 
trons "see" the quantum mechanical average value of a x , 
i.e. zero. Therefore the system behaves as if there was 
no spin present, i.e., one should substitute in Eq. © 
il — > (fi) = To- In the opposite regime, min[V,T] ^> A, 
the spin is in the mixed state and the electrons some- 
times "see" the spin in the state | f x ) and sometimes in 
the state | J. x ). Thus the current is the average of Ii x /i x - 
To find the output noise spectrum we perform the 
Laplace transform of Eqs. 1)25(1 

(J 2 )=Tr [[m u ,A} {s) p st ] 

+Tr [{[m u ,Z] {s) U(s)A(s)}p st ] 

-Tr [{{m u ,E] is) U(s)B(+0)} p st ] , (32) 

and, then, find the Fourier transform of the current- 
current correlator. This last step is done using Math- 
ematica as the expressions are quite extended. The 
Laplace transforms of the correlators a(t) and (3(t) con- 
tain, as usual, the real and the imaginary parts. It is pos- 
sible to show that the imaginary part of a gives rise to the 
renormalization of the system parameters, e.g., the Lamb 
shift of the qubit 's level splitting A, while the imaginary 
part of (3 is only important at very high frequencies and 
it ensures the causality of the meter's response functions 
(see Appendix O) . In what follows we neglect the Lamb 
shift as it is small compared to A. The full expression 
splits into three parts (I 2 ) = C\ + C-2 + C3 . 

The first term of Eq. (|32|l does not contain the evolu- 
tion operator U(s). It is, thus, expected to give a non- 
resonant contribution to the current-current correlator, 
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i.e., the pedestal (shot noise): 



Ci = 3o (s(w) +uj) 

s(uj + A) + s(w - A) + 2w 



-31- 



A (2A + s(w + A) - s(w - A)) 
2i(A) 



(33) 



For the symmetrized noise power Sj (uj) = Ci(w) + 
Ci{— to) we then obtain 



5f ot (^ = 0) = 2 5o s(0) + 2. 9lS (A)(l 



A 2 



(34) 



The two terms in Eq. (|34|l clearly correspond to the two 
terms of Eq. (pTT)l . This is a usual situation for the shot 
noise. Yet, the Fano factors for these two contributions 
are slightly different (at T = both are equal 1) . 

The last two terms of Eq. i1-'S2[) contain the evolution 
operator U(s) and, thus, are expected to produce reso- 
nant contributions. We, first, employ the Bloch-Redfield 
approximation, i.e., we substitute the self-energy £(s) in 
Eq. (|13() by Ebr- Then we obtain the two remaining 
contributions Ci and C3. The contribution C2 is 



C 2 (w) = 



(W) 2 rA 2 



{to 2 - A 2 ) 2 + 4r 2 cj 2 ' 
Ao(A) + wa(w) 



2Fs(A) 



(35) 



where 



5is(A) 



s(A) 



(36) 



is the qubit's dephasing rate, while 51 = I 
AV-y/gogi is the sensitivity of the meter. 

In the regime T <C A we obtain two peaks placed 
around to — ±A. This is how the qubit's damped coher- 
ent oscillations are reflected in the output noise. Inter- 
estingly, this contribution is symmetric, even though the 
symmetry should not have been expected in general. An 
asymmetric contribution would correspond to a change 
in the current-current susceptibility (finite frequency dif- 
ferential conductance) due to the presence of the spin. 
Such corrections have usually the Fano shape. As shown 
in Appendix [U] the lowest order Fano resonances vanish 
in our model due to the non-universal behavior of the 
QPC's response functions. 

The contribution C2 (Eq. I|35[ll is a product of a 
Lorentzian and a reduction factor in the brackets. The 
Lorentzian coincides with the one obtained in Refi .. The 
reduction factor simplifies for T = 0. Then, if V > A, 
it is given by (1 — A 2 /y 2 ), while for V < A it is equal 
to 0. In the last case the measuring device can not pro- 
vide enough energy to excite the qubit and, therefore, 
the qubit remains in the ground state and does not pro- 
duce any additional noise. The ratio between the peak's 



hight and the pedestal's hight is different for positive 
and negative frequencies. In the limit 171 <C go, T = 0, 
and A < V we obtain C 1 (±A) « g V(l ± A/V) and 
C 2 (±A) M 4 5o ^(l - A 2 /V 2 ) and, thus, 



C 2 {io = ±A) 



±A> - 4(1 4 } 



(37) 



For A — > V the ratio for the negative frequency peak 
reaches 8. In this limit, however, the peak's hight is zero. 
For symmetrized spectra the maximal possible ratio is 4 
(RefiiS). We see that the enhancement of the signal-to- 
noise ratio is due to the fact that the suppression of the 
qubit's contribution to the noise by factor (I — A 2 /V 2 ) 
is smaller than the suppression of the negative frequency 
background noise by factor (1 — A/V). An interesting 
question is what exactly is observed in the experiments. 
If the setup for the measurement of the noise would be 
absolutely passive, like the photon counters in the flu- 
orescence experiments, it could measure only what the 
system emits, i.e. the noise at negative frequencies^^. 
Moreover, if one is only interested in the excess noise, 
i.e., in the nonequilibrium addition to the noise power 
due to the finite transport voltage, then even an active 
detector may be useful. Namely, as shown in Refi^ 3 ., if 
the excess noise power is (almost) symmetric, it can be 
effectively measured by a finite temperature LCR filter. 
In our case the excess noise consists of the shot noise, 
C\ — Ci(V = 0), and the coherent peaks C2. While the 
second contribution is symmetric, the first one is only 
approximately symmetric in the limit g\ <C go- Thus, 
in this limit, the combination C2 + C\ — C\(V — 0) can 
be measured. The question of what can be measured in 
the regime g\ ~ go, when the excess noise is essentially 
asymmetric, will be considered elsewhere. 

Note also that far from the resonance, for ^> V, the 
reduction factor in Eq. I|35ll becomes negative, creating 
a very small negative contribution to the current-current 
correlator. This is an artifact of the Bloch-Redfield ap- 
proximation. The non-Markovian corrections are ex- 
pected to compensate this negative contribution so that 
the correlator is positive and vanishes at high negative 
frequencies. 

As the voltage increases so that T ~ A, the peaks given 
by Eq. 1135(1 start to overlap. Then they form a single peak 
around lu = which starts getting narrower. Finally, 
when r 3> A, the width of the peak scales as A 2 /r -C T. 
This is the strong measurement or Zenoi^ regime. The 
meter manages to almost localize the qubit in one of the 
eigenstates of the measured observable (<r x ). The rare 
transitions (flips) give rise to the "telegraph" noise peak 
around lu — 0. The stronger is the measurement (r) the 
longer is the average time between the flips (Zcno effect) 
and, thus, the narrower is the peak. In this regime the 
output is classical and the reduction factor plays no role 
(is equal to unity). 

Finally, for the last contribution C3 we obtain 



C 3 (w) 



r 3 



4r 2 
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-A A co 

FIG. 5: The output noise for the following parameters: V = 
1.5A, T = 0.1A, go = 0.01, gi = 0.0064. 

[a(A + w) + a(A - oj)] 

^(A) X 
(4VAs(A) - 2A 2 a(A) - s 2 (A) [a(A + u) + a(A - u)] 

-As(A) [a(A + w) - a(A - w)] ) . (38) 

This contribution corresponds to a small peak around 
w = 0. For T < A, T = 0, uj < A, and V > A we obtain 

4r 3 A 2 / A 2 \ 

while for the same conditions but V < A the contribution 
C3 vanishes. To understand the physical meaning of the 
peak at zero frequency we note, that due to the asymme- 
try of the correlators a and (3 the expectation values of 
the current corresponding to the two eigenstates of the 
spin's Hamiltonian | | z ) and | [ z ) are different (classically 
they would be equal as in both states (<7 X ) = 0). Indeed 
substituting into Eq. (|28[) the density matrices | Tz)(tz I 
or I I-,) (iz I instead of p Bt , that is forcing the steady state 
to be one of the eigenstates we obtain for the respective 
currents I — goV + g%(V =F a (A)). As the qubit is coher- 
ent (under-damped) the back-action noise causes random 
transitions between the eigenstates of the qubit 's Hamil- 
tonian. This translates into the "telegraph" noise of the 
current. The effect is governed by the ratio A/V and is 
small in the limit when F » A. In the quantum Zeno 
regime (r 3> A) the contribution C3 is always negligible 
as compared with C\ and C2 ■ 

In Fig. El we plot an example of the output noise (I 2 ). 



V. CONCLUSIONS 

We have calculated the output noise of the point con- 
tact used as a quantum detector for arbitrary voltage 
and temperature. In the regime eV ~ A and T -C A 
the output noise is essentially asymmetric. The qubit's 



oscillations produce two peaks at u> = ±A. The peaks 
have almost equal height and, therefore, the negative fre- 
quency peak is much higher relative to it's pedestal than 
the positive frequency one. The peak/pedestal ratio can 
reach 8. As the negative frequencies correspond to emis- 
sion, this could be observed by further passive detectors. 
We have also obtained a "telegraph noise" peak around 
lu = for a purely transverse coupling. This peak ap- 
pears due to the quantum asymmetry of the noise spec- 
tra. It means that the detector discriminates not only be- 
tween the eigenstates of the measured observable {a x in 
our case) but also between the states of different energy. 
The results of this paper are obtained for the simplest 
and somewhat artificial model of a quantum detector. In 
particular, in this model, the leading contribution of the 
spin to the output current correlator is of the peak type 
at all voltages and vanishes at V = 0. In general this 
should not be the case, as the coupling to an additional 
(discrete) degree of freedom usually changes the response 
functions of the continuum (Fano resonances). In our 
system, however, certain properties of the response func- 
tions of the meter (see Appendix Q prevent the Fano 
resonances from appearing in the leading order of the 
perturbation expansion. It would be interesting to per- 
form analogous calculations for more realistic detectors 
like SET's or QPC's with open channels. 

Recently, Bulaevskii, Hruska, and Ortiai^ studied the 
problem of a spin in a magnetic field interacting with tun- 
neling electrons with arbitrary spin polarization. They 
considered the case of low dissipation, r <C A, and the 
tunneling electrons were coupled to all the projections of 
the spin operator. 
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APPENDIX A: APPLICABILITY OF THE 
MACDONALD'S FORMULA 

For classical currents it is convenient to use the Mac- 
Donald's formula to calculate the noise power of the cur- 
rent. Recently this formula has been applied in Ref4& to 
calculate the output noise of a QPC used as a measuring 
device (the same system as in this paper). As only the 
limit F> A was considered, the output signal was clas- 
sical and the calculation using the MacDonald's formula 
was well justified. In this Appendix we clarify whether 
this approach is applicable when the output is quantum. 
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The MacDonald's formula reads 



Si(u)) = 2uj 



d(m 2 (t)) . 



dt 



sm(ujt)dt 



(Al) 



where a(t) = (m 2 (t)) is the dispersion of the integral of 
current m = L I(t')dt' . One starts counting the charge 
that have tunneled starting from t — 0. One also assumes 
that at t = the spin's density matrix is the stationary 
one (see Eq. (|26p) and that m = at t = 0. Then, since 
<j(t = 0) = 0, we obtain 



Si(uj) 



[cr(s = iuj + 0) + cr(s = -iw + 0)] . (A2) 



To obtain a(s) we apply the propagator U(s, Ak) (see 
Eq. I^nj) to the stationary density matrix and, then ap- 
ply twice the operator m = id/dAk (as we do it after 
the propagator, i.e. at the left-most end of the Keldysh 
contour, we do not distinguish between m u and m d ). As 
a result we obtain 



a(s) = -Tr 



0- 



dAk 2 



U(s,Ak)p s 



(A3) 



The derivative over AA: in Eq. 
using Eq. lfT3|): 



d 2 



dAk' 



-U = UT,"U + 2UH'UT,'U 



3J can be calculated 



(A4) 



is easy to obtain iS' = 
m u , E] = A - B (see Eq. (J23J) and 



From Eq. © 
[(ia/5Afc)E] A ^ 

(i 2 )E" = [(i<9/<9Afc) 2 E] Afc ^ = [m u , A] + /i.e.. After 
some algebra we conclude that the MacDonald's formula 
gives for the noise Si an expression very similar to the 
(symmetrized) Eq. I|32|l . However, while in the last line 
of Eq. (EU there is B(+0), the MacDonald's formula 
puts B(s) into that place. With this substitution we 
obtain for the C2 contribution 



(«) 2 rA 2 



v ; (uj 2 - A 2 ) 2 + AT 2 uj 2 V 
(A + oj)a(A -oj) + (A- uj)a(A + uj) + 2uua(u) 
Ws{A) 

(A5) 

In the classical limit A/V — > the result (|A5|) coincides 
with the one obtained in full quantum mechanical calcu- 
lation H35[). However the corrections (even when A/V is 
small) are not reproduced. Indeed at T = and A < V 
the reduction factor in Eq. (|33|) is (1 — A 2 /V^ 2 ) while 
Eq. gives (1 - A 2 /2V 2 ). 



APPENDIX B: BEYOND THE 
BLOCH-REDFIELD APPROXIMATION 

We have also calculated the output noise without using 
the Bloch-Redfield approximation, i.e., substituting the 



non-Mar kovian self energy E(s) into Eq. (|13fl . In the 
regime r« A we found the following correction around 
uj = ±A 



6C 2 (lu 
2(3o 



±A) 



4g ogi VA(Lu 2 - A 2 



(uj 2 - A 2 ) 2 + 8A 2 r 2 
gi )VA + gi a(2A)( s (A)TA) 



5 (A) 



(Bl) 



This feature has a "derivative" (Fano) shape and it makes 
the main peaks (Eq. IJ3SJ) a bit asymmetric. As far as 
the coupling constants (conductances) are concerned, the 
correction l|rJlfl seems to be of the same order as the terms 
Ci and C3. Thus, the question arises, what does the 
Bloch-Redfield approximation exactly mean. Analyzing 
this question deeper we note, that within this approxi- 
mation the self-energy E and, consequently, the evolution 
operator U factorize into two parts: the part describing 
the diagonal (in the eigen-basis of the qubit's Hamilto- 
nian) elements of the density matrix (two modes with 
eigenfrequencies around oj = 0) and the part describing 
the off-diagonal elements (two modes with eigenfrequen- 
cies around uj = ±A). The first part is responsible for 
the contribution C3, while the second part in responsi- 
ble for C-x. The non-Mar kovian corrections to E couple 
these two pairs of modes. These corrections are, however, 
proportional to the deviations from the eigenfrequencies, 
e.g., to (uj — A) . Thus, they vanish exactly at the eigen- 
frequencies. More rigorously, since the width of the res- 
onances is proportional to T ~ gx, the non-Markovian 
corrections carry an additional factor of go or g\ within 
the resonances. There, the Bloch-Redfield approximation 
is well justified and the corrections are of the higher or- 
der in go,gi- Outside the resonances the Bloch-Redfield 
approximation may be not justified. There, however, the 
main contribution is the shot noise term C\ which does 
not depend on U(s) and is not sensitive to the Bloch- 
Redfield approximation. In Appendix [21 we will show, 
that, indeed, the Fano shaped contribution i|Bl|) is of the 
higher order than those, that could, in principle, have 
appeared together with the main peaks (|35|1 . 



APPENDIX C: STANDARD KELDYSH 
CALCULATION WITH MAJORANA FERMIONS 

It is possible to obtain Eqs. H33fl and l|35|l using the 
standard Keldysh diagrammatic technique 46 . For the 
two-level (spin- 1/2) system it is convenient to employ the 
mixed Dirac-Majorana-fermion representation (see e.g. 
Re£% 



<r+ = Vxf 

a z = l- 2/t/ , 



(CI) 



where / is the Dirac fermion while rj z is the Majorana 
fermion (i] z = (g + gt), so that {r] z ,r) z } = 2 (g being 




G, 



a 





V 



./ 



FIG. 6: The graphical representation of the Green's func- 
tions G a (Eq. (EU) and Gp (Eq. JC4j). The loops are here 
to remind that each of these Green's functions is actually a 
combination of two electronic ones. In the approximations we 
use the electronic lines appear only in such combinations. 



FIG. 7: The lowest order non- vanishing diagrams for the self- 
energies and ~E V . The solid lines are electronic ones, the 
dashed line is the Green's function of the Dirac (/) fermion 
and the dotted line is the Majorana fermion's one. 



another Dirac fermion)). 

Our purpose is to calculate the correlator (I(t)I(t')) 
which can be presented as one of the components 
of the current-current Green's function Gi(t,t') — 
-i{T K I(t)I(t')). Namely 



(I(t)I(t')) = i[Gi} 21 = iG> 



(C2) 



In what follows we use the (Keldysh) notations ex- 
plained in Ref4£. The current operator is given by 
/ = iil(X — X^), while the tunneling Hamiltonian (the 
vertex of the perturbation theory) is ifr = 0(X + X'). 
It is, thus, convenient to introduce the two following 
Green's functions: 



G a = 



-t(T K t [X(t) - XHt)} ■ i [X(t>) - Xt (*')]) 
-i(T K [X(t) + Xi(t)] ■ [X(t') + X^t')]) , 

(C3) 



and 



Gp = -i{T K i [X(t)-X^(t)] ■ [X(t') +X\t')]) , (C4) 

where the subscripts a and (3 point to an obvious re- 
lation to the functions a and (3 introduced above. In- 
deed G> = —ia(t — t') and G< = —ia(t' — t), while 
G| = (3(t - t') and Gj = -0{t' - t). The two lines 
of Eq. I|C31 might in principle be different, for example, 
in the superconducting case. In our case, however they 
are equal. For these two Green's functions we use the 
graphical representation shown in Fig. 

Finally we introduce the fermionic Green's func- 
tions. For the Majorana fermions we define G^ = 
~i(TKr) z (t)rj z (t')) . It is easy to obtain the bare Green's 
functions G> = — i and G< = i. For the / fermions 
it is convenient to use the Bogolubov-Nambu represen- 
tation, i.e., * = (/, p) T and * t = (/*,/) and, then, 
G* = -i{T K *(t)&(t')). 

The bosonic functions G Q and Gp describe the reser- 
voirs and, therefore, are well approximated by their un- 
perturbed values. Making this approximation we neglect 



a possibility of the spin-induced correlations in the reser- 
voirs (leads), i.e., the Kondo effect. This is justified if 
max(T, V)>T K . 

The fermionic Green's function G^ and G^ describe 
the spin, which can be driven far out of equilibrium. 
Thus, we find these functions from the the kinetic 
(Dyson) equations. 



(C5) 



(C6) 



where S^, and are the / fermion's and 77 fermion's 
self-energies respectively. All the quantities in Eq. (|05jl 
are matrices 4 x 4 (in the Nambu and Keldysh spaces). 
The operator G^^ is given by 



r -i _ 



w-A 

uj + A 

ui-A 

uj + A 



(C7) 



while for G v q we obtain 



uj/2 
lu/2 



(C8) 



In Eqs. (|G7(I and (|G8|I we have neglected the infinitesi- 
mal terms responsible for, e.g., causality of the Green's 
functions. These are no longer needed when the finite 
self-energies are taken into account. 

For the self-energies we take the lowest non-vanishing 
order approximation shown in Fig. \7\ As usually, in the 
Keldysh space the self-energies are presented as 



5]\T/ = 



v-R pit 







yA 



y, A 







(C9) 



From Fig.Qit is easy to conclude that EjjJ = ilfX G> 
and = zT^AG^G^, where the Nambu matrix 



G> 
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A is defined A = 



Analogously we ob- 



tain 



iT 2 G> (1 -l)G> 



and E^- 



T i G< ( 1 -1 ) G% ( _M. To calculate the self- 



energies in the lowest order we can use the unperturbed 
retarded and advanced fermionic Green's functions. This 
is not so for the Keldysh component, which contains the 
information about the distribution function: G^i^iyi) = 

hgf/ v (w) (G$ /n (u) - G^j n {uj)\. As was pointed out in 

Refits the distribution functions hmj v {uj) are determined 
by the reservoirs even in the zeroth order. Thus they 
should be found self-consistently. After some algebra we 
find 



- = Tf\ (G* + h v (0)(GX - G£)) (CIO) 



and 



=T*\{{G%-G£) + h v {Q)GZ 



(Gil) 



In what follows we will only need the Green's function 
Gp. Therefore, instead of proceeding with the self- 
consistent determination of the functions hq,/ n , we note 
that h v (0) = just by symmetry (the self-consistent cal- 
culation gives the same). Thus we obtain 

lmE$(u) = -XT U , (C12) 

where T w = gis(u>) (see Eq. (|36|) where we have in- 
troduced T = Ta). Analogously, ImS^(w) = \T(w), 
The real parts of the retarded and advanced self-energies 
give the non-equilibrium generalization of the Lamb shift. 
Here we neglect it. For the Keldysh component we have 



S^(w) = — 2igiXuj 



(C13) 



Substituting the self-energy into the Dyson equation 
(|C5|) we obtain 



Gq, 


1 Sl-X 








/ —iv u «r w 


—2ig\ui 


2igiu> 






2ig\ui 


~2ig\io 


1 











V o o 







(C14) 



which is easy to invert. As a result we obtain 



G 



GK 



R/A 



- A ± iT u 



- A ± iT u 



{uj 2 - A 2 ) ± 2iuT u 



2ig\ui 



(w + A) 2 
2 - A 2 



uj - 



A 2 
-A) 2 



{uj 2 - A 2 ) 2 + 4w 2 r 2 



(C15) 



(C16) 



Now we are ready to calculate the current-current cor- 
relator \Vj2\ . The lowest order diagrams contributing to 
this correlator are shown in Fig. |H1 They give 

SG> (t - t') = T 2 G> (t - + iT ?n> (i - t')G> (t - , 

(C17) 



where 



n = -i{T K a x {t)a x {t')) 

= -i(T K (nzf + f f Vz\t)(Vzf + /V)( t ')> -(C18) 

One obtains Eq. (|C17I) summing all possible orientations 
of the lines in Fig. [5] 

In the Majorana representation the Green's function II 
is a "two-particle" Green's function (a bubble). To calcu- 
late it properly one has to take into account, e.g., the ver- 
tex corrections, which seems to be complicated. Instead 
we use here an identity, recently proven in Refs^S*£i, 
which reduces II to a single-fermion Green's function. 



M*K(*')> = ([/(*) + f (*)][/(*') + P (*')]) ■ (ci9) 



From Eq. (|C19(I it is easy to obtain the following rela- 
tions 



n> = (i i)g>(; 



and 



n = ! i ! )G<( \ 



(C20) 



(C21) 



Performing usual Keldysh manipulations we obtain 
-SigiuiA 2 



n^-ir 4 



{uj 2 - A 2 ) 2 + 4w 2 r 2 



n 



K 



-st^a 2 



{uj 2 



.2\2 



4w 2 r 2 , 



(C22) 



(C23) 



and, finally, 



1 



n> = -m K 

2 y 



Il R - IT 4 ) = 



-4igi(s(u}) + uj)A 2 



{uj 2 - a 2 ) 2 + 4^ 2 r 2 ■ 

(C24) 

We then obtain Eq. |j33J| from Eqs. (|GT7|l and lfU24|) in the 

regime T <C A. In this case we can approximate Il > by a 
sum of two delta functions and perform the convolution 
in the second term of Eq. (|G17ll . 

The diagrams giving the peaks at uj = ±A are shown 
in Fig. They are chosen out of many other second- 
order diagrams, since only in these diagrams the spin's 
line (the combined loop of the / and Majorana fermions), 
which gives the Green's function II, carries the external 
frequency uj. In all other diagrams the spin's lines par- 
ticipate in loops and are, thus, being integrated over the 
frequency. Then the resonant structure is washed out 
and one merely gets a second order correction to the shot 
noise (pedestal). Although in Fig. [5] we draw the loops 
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a) V 



T,® 




b) T, 




FIG. 8: The diagrams leading to Eq. . The /-line should 
be found from the Dyson (kinetic) equation as in Refi— . 



(bubbles) of / and r\ z lines, we do not actually calculate 
those but use instead Eqs. I|C22|) and (|C23|> . 

In it quite easy to calculate the first diagram (Fig.^Ji). 
It is given by T^G p {u})Ii{uj)[-Gp(u})\. In the other 
three diagrams the internal vertices, over which the in- 
tegration is performed, are not further connected. Let 
us, for example consider the diagram b). Acting, first, in 
the (11,12,21,22) Keldysh coordinates (see Refi&) we 
see that the left electronic loop of this diagram, after 
integration over the time of the "free" vertex, gives a 
"Keldysh vector" Xj , where j is the Keldysh index of the 
left external vertex. Thus the whole diagram b) can be 
presented as T^T^XjDjk (no summation over j), where D 
denotes the rest of the expression which can be treated 
as a usual Keldysh matrix (D in this case is given by 
II[— G/s] in the (R,A,K) coordinates). For Xj we obtain 
Xj = [-Gp(0)]u - [-Gp(0)]v = ~Gf(0). Thus X j is 
actually a (Keldysh) scalar and the diagram b) can be, 
finally, calculated as T 2 T 2 [-Gj (0)] • II(w) [-G (ui)] . Col- 
lecting the rest of the diagrams we obtain the following 
contribution to the current-current Green's function 



d) 




FIG. 9: The diagrams leading to Eq. 1351 . The solid lines 
are electronic ones, the dashed line is the Green's function of 
the Dirac (/) fermion and the dotted line in the Majorana 
fermion's one. 



Then, the first term in Eq. (|C27|I is the noise of the spin 
being amplified by the QPC. The interpretation of the 
second term in Eq. (|C27|) is less trivial. Now these are 
the internal correlations of the QPC being amplified by 
the combined system of the spin and the QPC. This con- 
tribution is negligible when uj <C V but is of the same 
order as the first one for oj ~ V. Using Eqs. (|C22(1 and 
(|C23I) we, finally obtain the following contribution to the 
current-current correlator 



SGj = TqT([Gp(u) - Gf(0) ■ 1] IIH [G£(0) • 1 - Gp(co)] . 

(C25) 

We calculate the Keldysh components of the Green's 
function Gp for u -C D, where D oc p^/R 1S the electronic 
bandwidth (the Fermi energy). As a result we obtain: 
G£(w) = nV[l + iO(u/D)], G$(u>) = -r)V[l - iO(u/D], 
and Gf(cj) = 2r]a(u). The factors 1 ± iO(u/D) are 
responsible for making the functions G R (i) and G^(t) 
causal. As we are interested in the low frequencies 
(u> <C D) we approximate those factors by 1. Then we 
obtain 



5Gi = 4goffi 



V a{uj) 




n R n K 
o n A 



—a(ui) 
V 



and, finally, 6G^ A = and 



(C26) 



SGf = 2SG> = 4.90.9! V 2 [U 



a(uj) 



(n R - n A ) 



(C27) 

The first term in Eq. I|C27I) is the standard contribution 
obtained in the high voltage limit, e.g., in Refs(i2*22. It 
can also be obtained by treating the QPC as a linear am- 
plifier, i.e., assuming the relation I(t) = I (t) + ^-a x (t)^. 



iSGf = 



(<5/) 2 r w A 2 



A 2 ) 2 + 4T 2 w 2 



uja(uj) 
Vs{u) 



(C28) 



The contribution (|C28(I coincides with C^w) 
(Eq. 123)) in the limit T < A. 

Note that the response functions G^(uj) and G A (w) 
vanish at V = 0. This non- universal property makes 
the spin's contribution to the equilibrium output current 
correlator to vanish. At any voltage and temperature, 
the contribution (|C28|) is of the "peak" type, and not of 
the Fano type. The Fano shaped resonances could have 
originated from the combination II R + II A = 2ReIT R . 
Indeed, for T <C A and \lu\ w A we obtain from Eq. 1C22|I 



IP 



IF 



4A 2 



s(A) {lo 2 



Alu 2 T 2 



(C29) 



This combination does not appear in Eq. (|C28fl due to 
another non-universal property of the response functions 
G R (tt>) and G A (w) which are purely real up to the fre- 
quencies of order D. Had this not been the case, the Fano 
type contribution would be similar to that of Eq. (|B1|I but 
of a lower order in go, g\ (oc 5o5i)- 
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